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We investigate the behaviour of a three-phase composite, structured as a hexagonal array of coated 
cylinders, when the sum of the dielectric constants of the core and shell equals zero. In such cases, 
the absorption of the composite is the same as the absorption of a periodic array of solid cylinders of 
core material and radius equal to the outer radius of the original coated cylinder. When the sum of 
the dielectric constants of the shell and matrix equals zero, the composite has the same absorption 
as a periodic array of solid cylinders of core material, and radius exceeding the shell radius. 



I. INTRODUCTION 

Previous studies of the effective dielectric constant of 
a periodic array of coated cylinders, in the quasistatic 
limit, have revealed some unexpected and interesting re- 
sults [1~3]. Specifically, we have used the method of Lord 
Rayleigh [4] to investigate the analytic properties of the 
effective dielectric constant of the array, as a function of 
the core dielectric constant (e c ) and the shell dielectric 
constant (s s ), while keeping the matrix dielectric con- 
stant (e m ) fixed. 

We have shown that when e c + e s — 0, the composite 
has exactly the same effective dielectric constant as a pe- 
riodic array of solid cylinders with dielectric constant e c 
and radius equal to the outer radius of the original coated 
cylinder. We have also shown that when e s + s m = 0, the 
composite has exactly the same effective dielectric con- 
stant as a periodic array of solid cylinders with dielectric 
constant e c , and radius exceeding the shell radius. In 
both cases, the core is magnified and it may be possible 
for a system with a vanishingly small concentration to ex- 
hibit properties like those of a concentrated system. This 
is what we have called a partially resonant system, since 
the system has a finite, rather than an infinite, response 
to the applied field. 

The same effect appears in the case of a random mix- 
ture of coated cylinders with nonlinear cores and linear 
shells [5]. Under conditions of partial resonance, the sys- 
tem behaves as a mixture of nonlinear solid cylinders of 
radii equal to, or larger than, those of the original coated 
cylinders. 

Here, we show how these results carry over to prob- 
lems of electromagnetic scattering by periodic arrays of 
coated cylinders, and elucidate the relation between par- 
tial resonances and anomalous absorptance of metallic 
wires [6] . From the many possible ways of packing cylin- 
ders in regular arrays in two dimensions we concentrate 
on the hexagonal lattice that is isotropic for fields ap- 



plied in the plane perpendicular to the cylinder axes (see 
Fig. 1). However, the method we use may be applied 
equally well to any periodic array of coated cylinders. 




FIG. 1. The Wigner-Seitz cell (red) with the central in- 
clusion (coated cylinder with r c — 0.03d, r s = 0.31d and 
d — 1) and the nearest neighbours (located along the 
dashed circle) for a hexagonal array. The centers of the 
nearest neighbours are given by the lattice vectors R p for 
p e {(-1,0), (-1,1), (0,-1), (0,1), (1,-1), (1,0)}. We also 
show the fundamental translation vectors ei = d(l,0) and 
e 2 = d(l/2,V3/2). 

Recent work by Pendry et al. [7] has heightened inter- 
est in the extraordinary and sometimes paradoxical be- 
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haviour which can arise in composite systems, in which 
one component has a dielectric constant (or a magnetic 
permeability) which takes approximately a real and neg- 
ative value. In such a composite we may observe strong 
focusing of light at a plane interface, or ultra-refraction 
[8], or a dilute composite masquerading as a concentrated 
one [1-3] . Again, the use of metal-coated spheres in pho- 
tonic crystals has resulted in robust band gaps, indepen- 
dent of sphere-packing geometries [9] . 



II. THE FIELD IDENTITY 

A. Field Expansions 

We consider the electromagnetic modes of a three- 
phase composite consisting of an array of coated cylin- 
ders of infinite length, each aligned parallel to the z 
axis, inserted into a matrix of dielectric constant e m and 
magnetic permeability /i m . For each cylinder, the core 
and shell are characterized by the radii r c and r s , di- 
electric constants e c and e s , and magnetic permeabilities 
li c and ii s , respectively. We also assume that the elec- 
tric and magnetic fields depend on time through a factor 
cxp (— iut). The core and shell of all cylinders, as well 
as the matrix are homogeneous so that, the fields satisfy 
the Maxwell equations: 
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Here, in the xy— plane, fl c denotes the union of all the 
domains occupied by cylinder cores, while Cl s denotes the 
union of all the domains occupied by the cylinder shells. 
Thus, f2 c 1J Sl s represents the region occupied by all the 
cylinders in the array and il m is the region between the 
cylinders (the matrix region). 

By eliminating E or H between (1) and (2), we obtain 
the Helmholtz equations for the components of the fields: 
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where n 2 = u 2 e^ is again a function of position. The axes 
of the cylinders in the array are parallel to the z— axis, 
and the fields are taken to depend on z through the fac- 
tor exp(i/3z), with /3 representing the propagation con- 
stant along the z direction, in all the regions O c , fi s and 



Sl m . Also, in cylindrical coordinates (r, 9, z), the field 
components E z and H z determine the other four field 
components through the relations [10]: 
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The field components E z and H z satisfy the Helmholtz 
equations: 
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where is the transverse Laplacian and 
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represents the transverse (in-plane) propagation con- 
stant. 

For the cylinder centered at the origin of coordinates, 
we represent the electric and magnetic fields E z and H z 
(denoted here by V), by series expansions in terms of 
cylindrical harmonics: 
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where J and Y represent the Bessel functions of the first 
and second kind. The three forms of the series expansions 
in (12) correspond to < r < r c (i.e., inside the core), 
r c < r < r s (i.e., inside the shell) and r > r s (i.e., in the 
matrix), respectively. Also, the superscripts c, s and m 
label the fields inside the cylinder core, cylinder shell, and 
in the matrix, respectively. Thus, we have k 2 — uj 2 e c fi c , 
k 2 = Lu 2 e s [i s , and k 2 m = u 2 e m ^ m - 



B. Boundary Conditions 

The boundary conditions express the continuity of the 
tangential components of the electric (E z , Eg) and mag- 
netic (H z and Hg) fields across the core and shell sur- 
faces. When the crystal momentum is perpendicular to 
the axes of the cylinders in the array, we have /3 = 0, 
and the problem can be reduced to solving two indepen- 
dent problems [11]: (i) for TM or s polarization (when 
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H z = and the transverse parts of E and H are gener- 
ated by Vj_l? z ), and (ii) for TE or p polarization (when 
E z = and V j_H z gives the transverse components of E 
and H). Thus, depending on the polarization, we obtain: 
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for TE polarization. 

In the case of TM polarization, for the cylinder cen- 
tered at the origin of coordinates, we use the series ex- 
pansions (12) for E z , and by means of the boundary con- 
ditions (13-16) we obtain relations of the form: 
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In (21), the prime indicates the derivative of the corre- 
sponding function and Zi = ^/Jm/sI (i = m, s, c) rep- 
resent the impedances of the matrix, shell and core, re- 
spectively. For TE polarization we obtain the relation 
between AY" and BY" by changing Z^ — > 1/Zi in (21). 

The structure of (21) allows us a generalization to mul- 
ticoated cylinders, and for this reason we have added the 
second column in the last matrix. However, this last ma- 
trix is multiplied by the vector L4f,0] T corresponding to 
the field expansion inside the core (i.e., the first form in 



(12)). Note that the first two matrices in the right side of 
(21) arise from the boundary conditions at the interface 
between the shell and matrix ((15) and (16)), while the 
next two matrices come from the boundary conditions at 
the core - shell interface ((13) and (14)). Consequently, 
in the case of a multicoated cylinder, we may infer the 
form of a general term representing the interface i (be- 
tween the i th and (i + l) th media) characterized by the 
radius n: 
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Here, we have counted the layers from the exterior to 
interior and thus, for n — 1 shells, the matrix is indexed 
by i = and the core (the innermost medium) by i = n. 
Now, Eq. (21) takes the form 
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We also obtain the expression for the TE polarization by 
substituting Zi — > 1/Zi. 



C. The Generalized Rayleigh Identity 

For each polarization (TM or TE), we may write a 
Generalized Rayleigh Identity of the form [12,13]: 
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with Mi from (22) and — oo < I < oo. The quantities 5 ; y 
are dynamic lattice sums, which may be evaluated using 
accelerated summation over the reciprocal array [12]: 
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Here, n is an arbitrary non-negative integer, £ is the 
length of an arbitrary vector inside the unit cell (shorter 
than the shortest line connecting two of its vertices), and 
A represents the area of the unit cell. The vector k 
is the in-plane component of the "crystal momentum" 
kfl = ko + /3z, with z the unit vector along the z— axis. 
Note that, in this paper we consider only the case (3 = 0, 
so that ks = ko. 

The quantities Q/j and 9h from (26) depend on the type 
of the array. In our case, the hexagonal array is defined 
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by the fundamental translation vectors ei = c?(l,0) and 
e 2 = d(l/2, V3/2), where d is the array constant (see 
Fig. 1). Hence, the centers of the cylinders are located 
at the points given by the radius vectors 



R p =_piei + p 2 e 2 , 



(27) 



with p denoting the pair of integers {pi,p 2 ) € Z 2 . The 
fundamental vectors of the reciprocal array are given by 
the relations 
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where ^4 = ||ei x ©2 1| and e 3 = (ei x e 2 )/A is a dimen- 
sionless unit vector, along the z axis. With the vectors 
(28) and (29) we form the reciprocal lattice vectors 



K h = /iiUi + h 2 u 2 , 



(30) 



where h — (hi,h 2 ) G Z 2 . Finally, we have = k + 
and Oh = arg(Qh). 

The eigenfrequencies (u>) of the Hclmholtz equation (5) 
are obtained from the zeros of the determinant of (25). 
In the coordinate system u versus ks, the photonic band 
diagrams are given by the trajectories of the eigenfre- 
quencies lo when ke follows the boundary of one of the 
irreducible regions in the first Brillouin zone (see Fig. 2). 

In the case of metallic media, k^ becomes complex. It 
is then simpler to use a summation method [14] different 
from (26), and a technique, based on Bloch's theorem, 
developed originally by McRae [15] in low energy electron 
diffraction, and applied recently to photonic crystals by 
Gralak et al. [8]. We consider the array as an infinite 
stack of gratings and calculate the scattering matrices 
for a single grating in the array [16]. Then, the complex 
crystal momentum k^ is found by solving an eigenvalue 
problem [17]. 



III. THE LONG WAVELENGTH LIMIT 

A. Partial Resonances 

The solution of the problem of electromagnetic scat- 
tering by a solid cylinder, at normal incidence and in the 
long wavelength limit [18], revealed the occurrence of a 
resonance in scattering when the relative dielectric con- 
stant of the cylinder was -1. Later, the existence of this 
resonance was proved experimentally in studies of radio 
reflection by meteor trails [19], and of plasma columns 
[20,21]. Also, Wait [22,23] has shown that in the long 
wavelength limit the scattering resonances given by the 
formula 



where the label 1 refers to the cylinder and 2 refers to 
the homogeneous medium outside the cylinder, are inde- 
pendent of the polarization and angle of incidence. Note 
that, in electrostatics the condition (ei + e 2 ) = de- 
fines the accumulation point for the poles and zeros of 
the effective dielectric constant of a two-phase compos- 
ite [24,25], with (/ii + /i 2 ) =0 corresponding to magne- 
tostatics. 

We applied the method developed by Wait [22] to solve 
the problem of scattering of a plane wave by a coated 
cylinder at oblique incidence. In this case, in the long 
wavelength limit we obtained the resonance condition 

(fi c + n s ) {/j, s + Hm) (e c + e s ) (e s + e m ) = 0, (32) 

which describes precisely the electrostatic partial reso- 
nances of a three-phase composite [1-3] , plus their mag- 
netostatic counterpart. The resonance condition (32) 
also suggests a generalization to a n-phase composite, 
of the form 
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(33) 



for the case of a multi-coated cylinder. Here, fi n and e„ 
denote the magnetic permeability and dielectric constant 
of the innermost phase (core), while /iq and £o represent 
the magnetic permeability and dielectric constant of the 
medium outside the composite. 



B. TE Polarization 

We concentrate now on non-magnetic materials for 
which /i m = fi s = [i c = no, so that e m = u^Sq, s s = n 2 s s 
and e c = n 2 e , where e is the dielectric constant of free 
space, and (i = m,s,c) represent the refractive in- 
dexes of the matrix, shell and core, respectively. In the 
long wavelength limit (k m — > 0), we use the expressions 
of the Bessel functions for small arguments [26] in the 
boundary conditions coefficients M; from (22), for TE 
polarization. For ! ^ we obtain 



M l =M_ l = — 



(34) 



where 
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while for I = 0, we obtain a completely different form 
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(pi + Hi) (ei + £2) = 0, 



(31) 



To relate the long wavelength limit of the dynamic 
problem with the corresponding problem in electrostat- 
ics we apply the same method as in Ref. [27]. Thus, 
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the boundary conditions (17-20) correspond to an elec- 
trostatic problem in which the inverse of the dielectric 
constants (e c — ► l/e c , e s — > l/e s and e m — > l/e m ) have 
to be considered. This will also change ji — ► —7;. Now, 
the boundary conditions (22) for Z ^ can be written in 
the form 



7T V k 



Dm 

- 1)! T , -| r . (37) 



In electrostatics, the corresponding relationship be- 
tween the coefficients Ai and which controls the re- 
sponse of a coated cylinder to an external field, has the 
form [1] 



(38) 



Note that, due to the symmetry properties of the elec- 
trostatic potential, I cannot take even values in the case 
of a square array, or values which are multiples of 3 in 
the case of a hexagonal array [28]. By comparing (37) 
with (38) we may deduce the relation between static and 
dynamic multipolc coefficients 
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for I 7^ 0. The dynamic field in the matrix (H 
be approximated by [27,29] 

Hi m \r) = A cxp [ik B V m {v)\ + 0(k 2 B ) , 
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where V m is the solution of the static problem. Note 
that the derivation of (40) does not involve the boundary 
conditions coefficient M from (36). 

In electrostatics, the partial resonances of a three- 
phase composite consisting of coated cylinders are de- 
fined by the equations 

e c + e s =0 (core shell resonance), (41) 

s s + £ m = (shell matrix resonance), (42) 



when (38) becomes 
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respectively. In the first case, the field inside the coated 
cylinder is exactly the same as would be found within a 
solid cylinder of radius r s and dielectric constant e c , while 
the potential outside the coated cylinder, in the matrix, 
is precisely the same as that outside the solid cylinder 
[1] . The second case corresponds to an effective dielectric 
constant of the array of coated cylinders, identical with 



that of an array of solid cylinders of radius r 2 s /r c > r s 
(the geometrical image of the core boundary with respect 
to the shell outer boundary), and dielectric constant e c . 
Now, the field external to the coated cylinder and beyond 
the radius r 2 / r c is the same as that external to the solid 
cylinder [1]. 

Since it is the relationship between the coefficients Ai 
and Bi which controls the response of a coated cylinder 
to an external field, equations (38) and (37) show that 
this response is determined by 7; in electrostatics as well 
as in the long wavelengths limit of electrodynamics. The 
limiting process is smooth and, therefore, we expect a 
resonant behaviour of the array of coated cylinders, for 
nonzero frequencies, when one of the conditions (41) or 
(42) is satisfied. 



C. TM Polarization 

The behaviour of the array of coated cylinders is com- 
pletely different in the case of TM polarization. Now, 
in the long wavelength limit, the boundary conditions 
coefficients Mi from (22) take the form 



Mi = M_i 

= -(— 

for I ^= 0, and 
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, (45) 
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Following the same method as in Ref. [27] we consider 
the central equation (I = 0) of the Rayleigh identity (25) 

[M + S% (k m , k )] B + Sl n {k m , k )B n = , (47) 

and employ the forms of lattice sums for small k m [27] 
4 
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(48) 
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where the superposed bar denotes complex conjugation, 
A is the area of the Wigner-Seitz cell, # = ar g (ko) and 
ot = k m /k . We also introduce the relations between 
dynamic (B n ) and static (B' n ) multipolc coefficients [27] 
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( km ) 
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(51) 
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for n > 0. Here, the static multipolc coefficients are, to 
leading order, independent of k m . Finally, the equation 
(47) becomes 
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(So 



e s ){r c /r s ) 2 + (e s - e m ) A(l - a 2 ) 



B' 



0. (53) 



From all the other equations (I 7^ 0) of (25) we obtain 
B[ = so that, in the long wavelengths limit, to lead- 
ing order, the only non-trivial solution of the Rayleigh 
identity (25) is B' ^ and B' n = (Vn ^ 0), if 



— fc^c ~T" fs&s ~\~ fm^r 



(54) 



with f c = irr 2 /A, f s = n(r 



!)/Aand/ m = 1-nrj/A, 



denoting the area fractions of the phases inside the unit 
cell. From (12) we derive the form of the electric field in 
the matrix: 



(s c - e s ){r c /r s ) 2 + (e s 



B' . 



(55) 



Equation (54) represents the linear mixing formula [30] 
for the effective dielectric constant of an array of coated 
cylinders subject to a uniform field applied along the 
cylinder axis. Note that there are no terms of the form 
e c + £ s or e s + e m in (53) or (54) to indicate a core-shell 
or shell-matrix partial resonance. Again, the limiting 
process is smooth and, consequently, we do not expect a 
resonant behaviour of the array of coated cylinders, for 
any frequency, in the case of TM polarization. 



IV. NUMERICAL RESULTS 



section we will denote by e the relative dielectric con- 



stant e/sq 



(the square of the refractive index) . With 



this notation, we use the experimental data for the real 
and imaginary parts of the refractive index of aluminium 
oxide in the region 0.3/im < A < 20.0/im [31], to eval- 
uate the relative dielectric constant of aluminium oxide 
(e^j 2 Q 3 = n^y q ). Also, the relative dielectric con- 
stant of the second phase is Ed = 1.876. 

We assume that the composite is structured as a hexag- 
onal array (array constant d=l/im) of coated cylinders 
(core radius r c = 0.03/im, cylinder radius (core+coating) 
r s = 0.31/J,m) in vacuum. The Wigner-Seitz cell is shown 
in Fig. 1, while the irreducible part of the first Brillouin 
zone (KFMK) is shown Fig. 2. For this region of the first 
Brillouin zone, the coordinates of the corners are 



\3d'dV3j 



r(o,o), m o, 



2tt 
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A. Partial Resonances 

In the case of an absorbing (metallic-type) core and 
a dielectric shell, the constant, the positive value of the 
shell dielectric constant (e s = n 2 ,) makes impossible a 
relation of the type n 2 , + — corresponding to a 
shell-matrix partial resonance. Hence, we expect only 
core-shell partial resonances, which are defined by the 
equation 



1 A1 2 3 



= 0. 



(56) 



In this section we discuss the behaviour of a three- 
phase composite consisting of aluminium oxide (AI2O3), 
whose refractive index depends on the frequency, and two 
pure dielectric phases: one with a refractive index rid — 
1.37, and vacuum (no = 1). 




FIG. 2. The irreducible part (green triangle) of the first 
Brillouin zone. 

Note that, all the defining equations for partial reso- 
nances and the boundary conditions coefficients, in the 
long wavelengths limit or in electrostatics, are homoge- 
neous functions of dielectric constants. Also, note that 
the phases are non-magnetic materials. Hence, in this 
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FIG. 3. Real (solid curve) and imaginary (dashed curve) 
part of AI2O3 relative dielectric constant versus wavelength. 
The horizontal red line shows — while the vertical magenta 
lines mark the two solutions Ai and A2, of (57). 

Note that, n 2 , is real and so we will solve the real part 
of (56): 



3?e 



4 A1 2 3 



= 0, 



(57) 



and then check the magnitude of the imaginary part. 

In the wavelengths domain 0.3/im < A < 20.0/im, 
Eq. (57) has two solutions 
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Ai w 12.861/im, 



(58) 



Ai w 12.151/im, 



(63) 



giving a value of 



giving a value of 



n Al 2 3 



10" 



0.329 i. 



when substituted in (56), and 

A2 ~ 17.333/im, 

giving a value of 

n 2 MQ +n 2 d ^l0- 4 + 32.72i, 



(59) 



(60) 



(61) 



when substituted in (56). They are marked by the ma- 
genta vertical lines in Fig. 3. Actually, the magnitude 
of the imaginary part of (56) shows that we may have 
a partial core-shell resonance at Ai, while the second 
root appears only as the result of an anomaly of n^y q 

around A2. 

In the case of dielectric core and absorbing (metallic- 
type) shell we expect the coated cylinders to exhibit both 
types of partial resonances: core-shell 



4 A1 2 3 



= 



and shell-matrix 



1 A1 2 3 



0. 



The wavelengths at which the first type of partial res- 
onances occurs are given by Eq. (57), while the second 
type, is identical with the classical resonance [24,25] ex- 
hibited by solid cylinders from AI2O3 in vacuum. 
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FIG. 4. Real (solid curve) and imaginary (dashed curve) 
part of AI2O3 relative dielectric constant versus wavelength. 
The horizontal red line shows — = — 1, while the vertical 
magenta lines mark the two solutions Ai, A2, of (62). 

Thus, for the shell-matrix resonances, in the range 
0.3/im < A < 20.0/im, the equation 



^ e ( n Al 2 3 +" 
has two solutions (see Fig. 4): 



= 0, 



(62) 



1 A1 2 : 



+ < 



10 



-15 



0.234 i. 



and 



A 2 w 17.344/zm, 



giving a value of 

n Ah0 3 ' 



10 



-13 



32.611 i. 



(64) 



(65) 



(66) 



These wavelengths are marked by magenta lines in Fig. 4. 
Again, the magnitude of the imaginary part of (64) and 
(66) shows that we may have a partial core-shell reso- 
nance at Ax, while the second root appears only as the 
result of the anomaly of n^y q around A 2 . 
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FIG. 5. TE polarization: Reflectance (red), transmittance 
(green) and absorptance (blue) for a stack of solid cylinders 
of refractive index n c = nj^ 2 Q 3 anc ^ iac ^ ms Tc = 0.03/im 
(dashed curves), and for a stack of coated cylinders (solid 
curves). In both cases the stack consists of 30 gratings 
of cylinders, and represents a "slice" from a 2-dimensional 
hexagonal array. The vertical magenta lines mark the two 
solutions Ai and A2, of (57). 



B. Absorbing Core and Dielectric Shell 

Here, we consider that the coated cylinders have an 
aluminium oxide core (n c = n Al 2 3 )' a dielectric coating 
(n s = rid), and are embedded in vacuum (n m = 1). The 
coating is a lossless, dielectric material of constant (i.e., 
independent of frequency) refractive index n s — l.37. This 
is essentially the same as the refractive index of MgF 2 for 
1/mi < A <2.5fim [31,32]. 



1. TE Polarization 

In Fig. 5 we show a comparison between two types 
of stacks of gratings: a stack of 30 gratings of solid 
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cylinders having a refractive index n c — np i i 2 Q 3 an d ra ~ 
dius r c = 0.03^m (dashed curves), and for a stack of 
30 gratings of coated cylinders defined at the beginning 
of Sec. IV B (solid curves). In both cases, the period 
of the gratings is d = 1/im. Actually, we compare the 
stack of coated cylinders with a stack in which the cylin- 
der coatings have been removed (but keeping al the other 
parameters unchanged), to show the influence of the shell 
on the optical properties of the stack (reflectance, trans- 
mittance and absorptance) , for normal incidence and TE 
polarization. 

We can see in Fig. 5 an enhancement in absorptance 
for the stack of coated cylinders at Ai, while the stack 
of solid cylinders exhibits an anomalous absorptance at 
A c ~ 12.151/im < Ai, which corresponds to the classical 
resonance of a two-phase composite n^y q + = 

[24,25] (in fact, the wavelength for this classical reso- 
nance, A c = Ai, is given by (62)) Thus, in the case of 
TE polarization there is a partial resonance of the type 
core-shell at Ai = 12.861/im. At this wavelength, the 
anomalous increase of the absorptance of coated cylin- 
ders shows a lossy (or metallic) type behaviour, which is 
a characteristic of aluminium oxide and is never exhibited 
by the dielectric shell (lossless media). Note that the core 
radius is ten times smaller than the shell radius. We can 
also see in Fig. 6 that, at Ai the absorptance of the stack 
of coated cylinders (blue solid curve) is approximatively 
equal with the absorptance of a stack of solid cylinders, 
made from AI2O3, and having the radius r s = 0.31/im 
(blue dashed curve). 



R, T , A 




FIG. 6. TE polarization: Reflectance (red), transmittance 
(green) and absorptance (blue) for a stack of solid cylinders 
of refractive index n c = n Al 2 03 an ^ ra dius r s — 0.31pm 
(dashed curves), and for a stack of coated cylinders (solid 
curves). In both cases the stack consists of 30 parallel gratings 
of cylinders, and represents a "slice" from a 2-dimensional 
hexagonal array. The vertical magenta lines mark the two 
solutions Ai and A2, of (57). 

Note the oscillations of the absorptance and transmit- 
tance, for the stack of solid cylinders of radius r s and 
refractive index n c (Fig. 6). In contrast with this be- 
haviour, the absorptance and the transmittance of the 



stack of coated cylinders are almost constant, except a 
small region around Ai. If we label by cc the reflectance 
(R C c), transmittance (T cc ) and absorptance (A cc ) for the 
stack of coated cylinders, and by sc the corresponding 
quantities (R sc , T sc , A sc ) for the stack of solid cylinders, 
then, in the neighbourhood of the core-shell partial res- 
onance, we have: 

Rcc ~ T sc ~ 0, T cc ~ -Rsc, A cc ~ A sc . (67) 

Equation (67) shows that the equivalence between the 
two types of stacks is actually a sort of dualitym the sense 
that they have similar absorptance, but the resonant sys- 
tem resembles a material which absorbs and transmits 
the incident radiation, while the stack of AI2O 3 cylinders 
absorbs and reflects incident radiation. 

In conclusion, at Ai, for normal incidence and TE po- 
larization, the stack of coated cylinders with a very small 
metallic-type core, has really the same absorptance as 
a stack of solid metallic cylinders of an enlarged radius 
r s = 0.31/im. It looks like the core and the shell are 
replaced by a new material which has the same absorp- 
tance as the core material, an effect identical with the 
core-shell partial resonance from electrostatics [1-3]. 



R, T , A 




FIG. 7. TM polarization: Reflectance (red), transmittance 
(green) and absorptance (blue) for a stack of solid cylinders of 
refractive index n c = "-AI2O3 anc ^ ra( hus r c (dashed curves), 
and for a stack of coated cylinders (solid curves). In both cases 
the stack consists of 30 gratings of cylinders, and represents 
a "slice" from a 2-dimensional hexagonal array. The vertical 
magenta lines mark the two solutions Ai and A2, of (57). 

2. TM Polarization 

There is no partial resonance for a stack of coated 
cylinders, for TM polarization. Fig. 7 shows no change 
of optical characteristics around Ai. The second solu- 
tion of (57) is A2 ~ 17.333/im, and at this wavelength 
the stack of 30 gratings of coated cylinders exhibits an 
anomalous absorptance (see Fig. 7), which corresponds 
to the large imaginary part of (61). The stack of solid 
cylinders having a refractive index n c — n^y q and ra- 
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dius r c = 0.03/im, exhibits an anomalous absorptance 
(see Fig. 7) at A2, corresponding to (66). 

In both cases, the anomalous high absorptance is not 
a resonance effect, and it is produced by an anomalous 
behaviour of the relative dielectric constant of AI2O3. At 
A2 ~ A2 the real part of the relative dielectric constant 
of AI2O3 is relatively small, while the imaginary part 
reaches a maximum value of ss 30 (see Fig. 3). Therefore 
the absorptance of the core of coated cylinders, or of the 
thin cylinders of core material only becomes large (see 
Fig. 7). At the same time, for long wavelengths and 
TM polarization, the effective dielectric constant of an 
array of cylinders is given by the linear mixing formula 
(54). Thus, for an array of coated cylinders the relative 
effective dielectric constant is 

e* = f c — + f s — + f m , (68) 

while for the array of solid cylinders of core material and 
radius r c we have 

S* = fc— + fm ■ (69) 

Hence, for a pure dielectric shell (e s ^constant), in both 
cases, the behaviour of e* is determined by e c and this 
explains the common maximum in absorptance shown in 
Fig. 7, at the same position as the maximum in the imag- 
inary part of the relative dielectric constant of AI2O3 (see 
Fig. 3). Also, around the wavelength associated with this 
maximum, the real part of the relative dielectric constant 
of AI2O3 exhibits a sudden jump from negative to pos- 
itive values, which makes almost equal the roots of the 
equations (57) and (62). Actually, the anomalous ab- 
sorptance exhibited by stacks of thin, solid cylinders has 
been described in Ref. [6]. 

C. Dielectric Core and Absorbing Shell 

Here, we consider the complementary structure with 
coated cylinders, having a dielectric core (n c = rid) and 
an aluminium oxide coating (n s = n^y O3)' emD edded in 
vacuum (n m = 1). Now, the core is a lossless, dielectric 
material of constant e c . In this case, we expect both 
types of partial resonances to occur, and in both cases the 
coated cylinder has to exhibit the characteristics of a solid 
cylinder made from the core material. Absorptance is a 
characteristic feature of the metallic phase. Therefore, 
we expect a decrease of the total absorptance for the 
stack of gratings, around Ai and Ai. 

1. TE Polarization 

For TE polarization the reflectance, transmittance and 
absorptance of the stack of coated cylinders with dielec- 
tric core and metallic shell, are almost the same as in the 



case of a stack of solid cylinders of radius r s , made from 
aluminium oxide (see Fig. 8). Thus, in the neighbour- 
hood of the shell-matrix resonance at Ai w 12.151/im the 
two stacks have the same optical characteristics. In fact, 
Eq.(44) shows that the coated cylinders must have opti- 
cal characteristics similar to those of solid cylinders made 
from core material and of an extended radius a = r^/r c . 
In our case this means equivalent dielectric cylinders of 
radius a — 3.2/im ^>> d/2 — 0.5/im, i.e. we have a 
stack of intersecting cylinders forming a homogeneous, 
dielectric slab, and we cannot apply the formulation from 
Sec. II. In the neighbourhood of the core-shell reso- 
nance at Ai w 12.861/im, the reflectance of the stack 
of coated cylinders is diminished, while the absorptance 
is increased, compared with a stack of aluminium oxide 
cylinders (Fig. 8). At A2 « A2 the two types of stacks ex- 
hibit identical optical characteristics, and the local max- 
imum in absorptance appears due to the anomalous be- 
haviour of £^20 aroun d this wavelength (see Fig. 4). 
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FIG. 8. TE polarization: Reflectance (red), transmittance 
(green) and absorptance (blue) for a stack of solid cylinders 
of refractive index n c — ^Ai 2 o 3 an< ^ racuus r ' s ~ ^-31 (dashed 
curves), and for a stack of coated cylinders with r c = 0.03 
and r s = 0.31 (solid curves). In both cases the stack consists 
of 30 gratings of cylinders in a hexagonal array. The vertical 
magenta lines mark the solutions Ai and A2 of (57), and Ai 
and A 2 of (62). 

To avoid such situations when the resonant system 
cannot be modelled by the formulation from Sec. II, we 
consider a stack of coated cylinders with a larger core 
r c = 0.21. In this CELS6, cL shell-matrix partial resonance 
will lead us to an equivalent dielectric cylinder, having a 
radius a = 0.46/im < d/2. Now, in the neighbourhood 
of the shell-matrix resonance at Ai, the stack of coated 
cylinders exhibits a drop in absorptance, an increased 
transmittance and a very small reflectance, while a stack 
of aluminium oxide cylinders having the radius r s has a 
zero transmittance and equal reflectance and absorptance 
(see Fig. 9). The equivalent stack, formed by dielectric 
cylinders of radius a, has a zero absorptance by definition. 
For long wavelengths (A > d) it has a low reflectance and 
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a high transmittance, therefore, the only indication of 
the existence of a shell-matrix partial resonance is the 
nonzero transmittance exhibited by the coated cylinders. 




A[jum] 



FIG. 9. TE polarization: Reflectance (red), transmittance 
(green) and absorptance (blue) for a stack of solid cylinders 
of refractive index n c = nj,^ q 3 and radius r s = 0.31 (dashed 
curves), and for a stack of coated cylinders with r c = 0.21 
and r s — 0.31 (solid curves). In both cases the stack consists 
of 30 gratings of cylinders in a hexagonal array. The vertical 
magenta lines mark the solutions Ai and A2 of (57), and Ai 
and A 2 of (62). 
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FIG. 10. TM polarization: Reflectance (red), transmit- 
tance (green) and absorptance (blue) for a stack of solid cylin- 
ders of refractive index n c = n A\ 2 Q 3 and radius r a = 0.31 
(dashed curves), and for a stack of coated cylinders with 
r c — 0.03 and r s = 0.31 (solid curves). In both cases the 
stack consists of 30 gratings of cylinders in a hexagonal array. 
The vertical magenta lines mark the solutions Ai and A2 of 
(57), and Ai and A 2 of (62). 

In the neighbourhood of the core-shell resonance (at 
Ai), Fig. 9 shows a reduced absorptance for the stack of 
coated cylinders. Again, the equivalent stack, formed by 
dielectric cylinders of radius r s , predicted by (43), has 
a zero absorptance by definition, and so the comparison 
with a stack of aluminium oxide cylinders of radius r s , 
seems to be more natural. In this case, it is interesting 



to remark the appearance of the duality mentioned in 
Sec. IV B 1. The two stacks have the same absorptance 
and the reflectance of solid cylinders is equal with the 
transmittance of coated cylinders 



R, 



T 

1 c, 



T 



A,. 



(70) 



2. TM Polarization 

As in the case of TE polarization, when the core radius 
of the coated cylinders is small (so that a — r 2 s /r c > d/2), 
for TE polarization the reflectance, transmittance and 
absorptance of the stack of coated cylinders with dielec- 
tric core and metallic-type shell, are almost the same 
as in the case of a stack of solid cylinders of radius r s , 
made from aluminium oxide (see Fig. 10). In fact, there 
is no resonant behaviour in the whole range 10/im < A < 
20^m. 

In the case of a larger core of coated cylinders, the ef- 
fect of shell-matrix partial resonance is small for TM po- 
larization. Around A = Ai the absorptance of the coated 
cylinders is almost equal with their transmittance, and 
the reflectance is close to zero. The curves displayed in 
Fig. 11 show no sudden change around Ai (shell-matrix 
partial resonance) and Ai (core-shell partial resonance), 
but it is interesting to remark the duality (70) at Ai. 
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FIG. 11. TM polarization: Reflectance (red), transmit- 
tance (green) and absorptance (blue) for a stack of solid cylin- 
ders of refractive index n c = nj^ 2 3 an d radius r s — 0.31 
(dashed curves), and for a stack of coated cylinders with 
r c — 0.21 and r s = 0.31 (solid curves). In both cases the 
stack consists of 30 gratings of cylinders in a hexagonal array. 
The vertical magenta lines mark the solutions Ai and A2 of 
(57), and Ai and A 2 of (62). 



V. CONCLUSIONS 

We have analyzed the partial resonances of a three- 
phase composite, structured as a hexagonal array of 
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coated cylinders, at long wavelengths. For such struc- 
tures there are two types of partial resonances: shell- 
matrix and core-shell. We have shown that the shell- 
matrix resonance has a very small effect. In contrast, the 
core-shell resonance appears for TE polarization in both 
cases of coated cylinders metal-dielectric-vacuum and 
dielectric-metal- vacuum. Note that, for such resonances 
the stack of coated cylinders has a nonzero absorptance, 
so we cannot compare it with a stack of dielectric cylin- 
ders (in the case of dielectric-metal-vacuum composite) . 
However, the comparison with a stack of solid, metal- 
lic cylinders of radius r s shows equal absorptance and a 
rcflcctance-transmittance duality, for the core-shell res- 
onance. 

An interesting behaviour is exhibited by the composite 
in the vicinity of shell-matrix (Ai) and core-shell (Ai) 
resonances (see Figs. 6, 7, 9 and 11). In these regions, 
the composite has a very low reflectance (R cc w 0), while 
the transmittance and absorptance take similar values. 

Generally, a band gap is defined by R=l, and T=0, 
so that sudden changes in optical characteristics of an 
array of coated cylinders will change the photonic band 
diagram. In our calculations the partial resonances oc- 
cur at long wavelengths (i.e., in the region of the acoustic 
band) and we expect the photonic band diagram of the 
composite to show horizontal lines intersecting the acous- 
tic band. These ultrafast changes in the behaviour of a 
composite, due to the dependence of the refractive in- 
dexes on the wavelength, can provide the basis of future 
optical switches in photonic integrated circuits. 
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